Charge dynamics in an ideal cuprate C d2- x N a x Cu02Cl2- optical conductivity from 

Yang-Rice-Zhang ansatz 
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We theoretically investigate charge dynamics in weakly coupled CuOi planes of the cuprate 
Ca2~ x Na x Cu02Ch (CNCOC) using Kubo formula for optical conductivity in the underdoped 
regime. The spectral function needed in Kubo formula is obtained from an analytical form of 
electron Green's function proposed (ansatz) by Yang-Rice-Zhang (YRZ) for the underdoped cuprates 
based on their previous renormalized mean field theory and on the investigations of weakly coupled 
Hubbard ladders. Although to an unaided eye the results of the numerical calculation look very 
similar to that found experimentally in [K. Waku etal. 2004[l|]] but a careful examination with 
extended Drude formalism shows that YRZ ansatz for the calculation of optical conductivity is not 
sufficient to understand the charge dynamics in CuOi planes of the cuprate CNCOC. More physics 
is needed especially electromagnetic response from bound charges. 



PACS numbers: 74.72.-h; 74.25. Gz; 74.72.Kf; 78.20.Bh 



I. INTRODUCTION 

Physical investigations at very small length scale (sub- 
atomic, atomic etc..) are usually indirect. To expose 
or uncover the underlying physical reality one has to 
relay on indirect experimental probes. In simpler situ- 
ations few experimental probes can give one a consis- 
tent physical picture (for example spectral lines of the 
atomic spectra suggested the discrete energy level struc- 
ture of atoms). In more complex situations (as for ex- 
ample in high temperature superconductors) one has to 
deal with several experimental probes. Consistent mi- 
croscopic physical picture of a material can only be con- 
structed if results of several experimental probes are con- 
sistent with each other and comply with one universal 
picture. And that universal picture is the major de- 
riving force, and at the same time it also puts heavy 
demands on researchers of understanding an arsenal of 
complex experimental probes (Fig. [1]). In complex sit- 
uations, understanding reached based solely on few ex- 
perimental results may not be consistent with other ex- 
perimental results. One has to be careful. The situa- 
tion is analogous to the famous story of blind men and 
an elephant. Each blind man "feels" different part of 
the elephant and concludes that the elephant is like a 
snake (for a man who touches the trunk) or like a pillar 
(for a man who feels the legs ) etc. Only if they pos- 
itively council (without fighting with each other) they 
may reach on the "universal" picture of the elephant. 
At present, for the cuprate 
high temperature supercon- 
ductivity problem we do not 
have a consistent universal 
picture for the whole phase 
diagram consistent with all 
experimental probes (in any 
way comparable to the de- 
scription of ordinary superconductivity by Bardccn, 
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FIG. 1: An arsenal of complex experimental probes. 



Cooper, and Schricffcr). But there are marvelous 
attempts Cuprates are complex materials with many 
anomalous properties-not only unconventional supercon- 
ductivity but also unconventional normal state response. 
Some of the anomalous properties are as follows. D.C. 
resistivity of hole doped cuprates just above the optimal 
doping is linear in temperature T (reminder: for a good 
metal obeying Fermi Liquid Theory (FLT), p ~ T 2 , and 
p ~ T if phonons also contribute at T (T << Tu e bye))- 
Away from optimal doping resistivity shows a very com- 
plex behaviour and electron doped compounds do 
not show this linear in temperature behaviour of resis- 
tivity! Hall coefficient is temperature dependent (for 
a good FLT metal it should be temperature indepen- 
dent). For cuprates, Drude scattering rate turns out 
to be frequency and temperature dependent (for a good 
FLT metal it is constant (again if phonons do not con- 
tribute in the temperature regime of interest)). In the 
"normal" state above optimal doping real optical conduc- 
tivity shows ~ — behaviour (in a good FLT metal it is 
~ -js Drude behaviour) . NMR relaxation rate shows sub- 
stantial deviations from linear in temperature behaviour 
(T- linear Korringa type). High T c cuprates has smaller 
coherence length (roughly the size of the cooper pair) as 
compared to conventional (in accord with BCS theory) 
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superconductors. The cooper pair size turns out to be 
smaller than mean carrier-carrier separation, thus mean- 
field type approximations in cuprate problem are ques- 
tionable. There are many other anomalous properties 
without proper understanding^. Cuprates show variety 
of phases with temperature and doping and one of the 
elusive phases is the pseudogap phase (Fig. [2]). There 
are many views on pseudogap phase [H, but recently, a 
broad picture at a phenomcnological level of the pseu- 
dogap phase is proposed for the cuprate problem. It is 
contained in Yang- Rice-Zhang (YRZ) ansatz This 
is based on several inputs @. 

The key element of YRZ's phenomenological theory is 
an ansatz for the single electron propagator. This ansatz 
is the outcome of series of investigations of the proposers 
over several decades and is primarily based on their renor- 
malized mean field theoryQ and on the investigations 
of weakly coupled Hubbard ladders Q. YRZ ansatz has 
been applied successfully to understand several anoma- 
lous properties of the pseudogap phase @. This has been 
applied with success to the interpretation of ARPES (an- 
gle resolved photoemission) experiments on CNCOC [9| 
and good agreement is seen between the calculated hole 
Fermi pocket and experimental data[f|. YRZ ansatz also 
correctly reproduce particle-hope asymmetry seen in ex- 
periments of Yang etalJTO]. This has also been applied 
to AIPES (angle integrated photoemission) and quali- 
tatively reproduce the key features of the spectra fa [Tl|. 
For STM (scanning tunneling microscopy) spectra of con- 
stant quasi-particle energy contours on BSCCO[lH one 
sees a qualitative agreement with YRZ, although quanti- 
tative fits with experimental spectra are not claimed]^]. 
Raman Spectra is quite useful because it probes both the 
nodal (B 2g ) and anti-nodal (Bi 9 ) charge dynamics (l3|. 
Valenzuela and Bascones[l4j found that YRZ qualita- 
tively reproduce features seen in the spectra(l5j (in per- 
ticluar they deduce two energy scales, nodal and anti- 
nodal, with opposite dependence on doping, nodal scale 
decreases with underdoping while the anti-nodal one in- 
creases). YRZ has been applied to several other exper- 
imental results see for detail @. In regard to the mi- 
croscopic picture, one should clearly distin guis h between 
YRZ ansatz and preformed pair picture [l6| in which 
particle-hole symmetry is maintained. 

Our interest here is in the optical spectra and in the be- 
haviour of optical conductivity in the pseudogap phase- 
its variation with doping and temperature. Optical con- 
ductivity spectra of the cuprates is complex in compari- 
son to that of conventional superconductors. In conven- 
tional superconductors the gap is isotropic in momentum 
and one gets a clear signature of the gap in the spec- 
trum (In fact early pioneering experiments by Tinkham 
etal[12l gave support to the energy gap model of super- 
conductivity on which Bardeen highly relied and later 
on a very successful theory of electromagnetic response 
was put forward by Mattis-Bardeen [l8j ) . In Cuprates 
gap(s) are anisotropic and optical conductivity spectra 
becomes intricate. Recently, experimentally observed op- 
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FIG. 2: Phase diagram of cuprates. 

tical spectra of cuprates has been studied using the YRZ 
theory by Carbotte, Nicol and colleagues [l9j . Their in- 
vestigations support the view that YRZ ansatz qualita- 
tively reproduce low energy behaviour of optical conduc- 
tivity. 

Here, in the present investigation, we re- visit this prob- 
lem. Our results do not bring any good news for the 
applicability of YRZ ansatz to the low energy optical 
response of cuprates. We concentrated on a specific 
compound Ca 2 _ x Na x Cu0 2 Cl 2 (CNCOC). The results 
of our numerical calculation appears to be in qualitative 
agreement with what has been found experimentally inQ 
(to an unaided eye the results look very similar to that 
found experimentally). But a careful examination with 
extended Drudc formalism shows that YRZ ansatz for 
the calculation of optical conductivity is not sufficient 
to understand the charge dynamics in Cu0 2 planes of 
CNCOC. It seems that more physics is needed to fully 
understand the optical response. 

Rest of the paper is organized as follows. In the next 
subsection (subsection A) essential points of YRZ ansatz 
and the cuprates are given. In subsection (B) a brief in- 
troduction to optical conductivity and Kubo formalism 
is given. In section II, CNCOC system and experimen- 
tal results are summarized. The calculation of optical 
conductivity using YRZ ansatz and comparison with ex- 
periment is given in section III. In section IV extended 
Drude model analysis is presented to point out the in- 
consistencies. We end with brief conclusion in section 
V. 



A. Essential points of YRZ ansatz and cuprates 

The starting point of YRZ ansatz is the reasonably well 
understood regimes of extreme underdoping and over- 
doping (Fig. 2). At overdoping one observes full Fermi 
surface and it disappears at zero doping. As one shifts 
from overdoping to underdoping in the phase diagram 
Fermi surface evolves from full Fermi surface to only 
disconnected arcs (in specific directions in fc-spacc) at 
underdoping to no Fermi surface at zero doping. The 
key issue is: how can one explain this doping depen- 
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dent evolution of the Fermi surface using a microscopic 
model? The major problem is the intermediate regime 
and YRZ ansatz is an attempt to fill this gap. As is 
well known, undoped cuprates are Mott insulators (more 
precisely charge transfer type) and it was pointed out 
by Phil Anderson very early on that the operational el- 
ements in the cuprate superconductivity are the CuOi- 
planes. Cu atoms in CuC>2 planes are in d 9 configura- 
tion (with one hole in the higher energy anti-bounding 

Zd x i_ y iCu 2Px(y) orbital lying in CuO^ plane). The 

holes localize (immobile) on atomic sites due strong on- 
site coulomb repulsion. With further hole doping, new 
holes are created in CuOi planes. These new holes will 
not be in Cu d-orbitals (because of strong Coulomb re- 
pulsion) and tend to be in oxygen p orbitals. It was 
shown by Zhang and Rice [2(| that if they form singlet 
pairs with the holes in Cu atoms then they will have 
lower energy. These singlet pairs are now called Zhang- 
Rice (ZR) singlets. We are interested in the transport 
properties of these ZR singlets. This also reduces the 
cuprate problem from three band to one band problem 
(although there are debates in the literature[2l[). The 
idea is to reduce the problem to its bare essentials. This 
motivates the famous t — J modclQ with no double occu- 
pancy at the same site (the Hilbert space of the problem 
will not have any configuration in which any site is dou- 
bly occupied-the projected Hilbert space) . One can work 
in unprojected Hilbert space by using the ideas of Renor- 
malized Mean Field Theory (RMFT)@, H- This is the 
first component of YRZ ansatz. The second component 
is based on the analytical studies by Konik, Rice, Tsvclik, 
and Ludwig 22] of 2-lcg Hubbard ladders and a collection 
of weakly coupled Hubbard ladders. For weakly coupled 
2-lcg Hubbard ladders the coherent part of Green's func- 
tion turns out to be (treating the inter-ladder coupling 
in RPA): 



G(k x , k y ,uj) 
G°(fc x ,w) 



1 



Gq 1 (k x , U>) — ty(ky) 
1 



(1) 



Where t y is the transverse inter-ladder hopping. e{k y ) 
and A are the bare band dispersion and quasi-particle 
gap respectively. 



Cuprates 



Many weakly coupled 2— leg Hubbard ladders 

FIG. 3: YRZ conjucture: equivalence of many weakly coupled 
2-leg Hubbard ladders and CuOi planes in cuprates. 

Based on the above two main elements, here is the key 
postulate: continuous crossover from the weak to strong 




interaction limit. The coherent part of YRZ Green's 
function for the pseudogap state of cuprates is postulated 
to be: 



G YRZ (k^) 
S Pff (fc, w) 



9t(x) 



huj — e(fc) — Y, pg {k,uS) 
|A PS (fc)| 2 



tun + eo(k) 



(2) 



This is essentially based on intuition (Fig. [3]). The 
dressed band dispersion e(k) contains the renormaliza- 
tion factors of RMFT. Details of bare band and renor- 
malized band dispersions are given in the next section. 



B. Optical Conductivity 

Optical conductivity is the linear response function of 
an external A. C. electric field. In a rough classical pic- 
ture the charge carriers oscillate back-and-forth under 
the influence of an external A.C. electric field (neglect- 
ing the weaker (-th) magnitude magnetic field effects). 
During this back-and-forth motion they scatter by im- 
purities and phonons causing the dissipation of energy 
(Joule heating). In the linear response regime induced 
current density J(r, uS) is related to the applied electric 
field E(r',w) by 

J(r,w) = y*cr(r,r';w)E(r',a;)dr'. 

If we are not in the anomalous skin effect regime 
(i.e., electric field does not vary substantially on the 
length scale of the mean-free-path of charge carriers) then 
er(r, r'; ui) = er(r; u))5(r — r') and 

J(r, w) = cr(r; w)E(r, lo). 

And further if the material is homogeneous: 

J(r, lo) = <x(w)E(r, u>). 

The optical conductivity ("Optical", if w is in opti- 
cal frequency range) cr(w) can be calculated by using 
quantum mechanical expression for current density J = 
a^iW^VV' - ipVtp*) and it can be shown that 



cr ( w ' T ) = -TTF v o,x( k *> k 



(3) 



dyfjL + Ml A(k x ,k y ,y)A(k x , ky,u> + y) 



This is called the Kubo formula (see for details [2J]). Here 
V is the volume of the sample, vo, x {k x , k y ) = j- de (^- fc «) 
is the Fermi velocity, f(u>) — / 3( fi ^_ Mp)+1 is Fermi-Dirac 
distribution function (/? = jjpf j P"v ^ s the chemical poten- 
tial). Spectral function A(k x , k y , uj) is given by the usual 
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expression A(k x , k y ,u>) — — ^ImG YRZ (k x , k y , oj + i0 + ). 
The coherent part of YRZ Green's function is postulated, 
as given in equation ([2]) , thus the spectral function takes 
the form: 



A(k x ,k y ,uj) = 



7 



)2 +7 2 



(4) 



Y*(k x , k v ,uj) is given in equation©, and broadening 7 is 
introduced phenomenologically[19|. This gives finite life 
time of excitations due various scattering mechanisms 
during transport. Rcnormalizcd dispersion e{k x ,k y ,oS) 
and bare band dispersion eo(k x , k y , u) in the YRZ 
theoryd Gj| are 

e(k x ,k y ) = — 2t(x)(cos(k x a) + cos(k y a)) 

— 4t'(x) cos(k x a) cos(kya) 

— 2t" (x)(cos(2k x a) + cos(2k y a)) — fi p , (5) 

and eo(k x ,k y ,Lu) = —2t(x)(cos(k x a) +cosffc y a)). Here a 
is the lattice constant. In YRZ theory [l, Q various band 
parameters are 

t(x) = g t (x)t + p{g s (x), £ = 0.338, J = |, 



t\x) 
t"(x) 



g t (x)t' , t' = -0.3to, 9t(x) = 



9t(x)t 



*o 



0.2to, g s {x) 



2x 

1 + x 
4 

(1 + z) 2 ' 



(6) 



as calculated from the band structure of Ca 2 Cu0 2 Cl 2 )Q, 
[24|. gt(x) and g s (x) are renormalization factors in 
YRZ theory. And the pseudogap A pg (k x ,k y ) = 
^A° g (x)(cos(k x a) - cos(k y a)) with A° g (x) = 0.6* (1 - 
x/0.2) in eV. Usually t ~ 400 meV in Cuprates. Be- 
fore we present the numerical calculation it is important 
to consider the experimental results of the system under 
consideration. 



II. THE Ca 2 - x Na x Cu0 2 Cl2 SYSTEM AND 
EXPERIMENTAL RESULTS 

The Ca 2 - x Na x Cu0 2 Cl 2 (CNCOC) system is a near 
ideal cuprate with single Cu0 2 plane per unit cell and 
the coupling of Cu0 2 -p\&nes in CNCOC is expected to be 
weaker than that of other cuprates like La 2 - x Sr x Cu0^ 1 
YBa 2 CusOT, Bi 2 Sr 2 CaCu 2 Os etc. This is due to the 
fact that in CNCOC the apical atoms are Chlorine ions 
(instead of Oxygen as in other cuprates) and Chlorine 
planes have more ionic character. This exert pulling of 
the d-electron on Cu ion towards the apical axis and 
thereby leads to very pointed octahedron (more positive 
charge on the Cu ion pull the planner Oxygen atoms to- 
wards itself and thus pushing apical atoms further out 
and creates a pointed octahedron). Thus the c-axis dis- 
tance in CNCOC is about 15 A, larger as compared to 
that in other cuprates (~ 10 A). Also CNCOC has no 
orthorhombic distortion from tetragonal structure as it 



is cooled through the pseudogap boundary. Due to these 
qualities CNCOC can be regarded as a better cuprate 
regarding charge dynamics in Cu0 2 planes. 

Charge dynamics of CNCOC has been measured exper- 
imentally in a beautiful piece of work by Waku etal.[l|. 
They measure D.C. resistivity and A.C. optical con- 
ductivity of single crystals of CNCOC grown by flux 
method [l|. They also measure temperature dependence 
of low energy (hco < 1 eV) optical conductivity at var- 
ious doping levels (Mott gap appears at higher energy 
thj ~ 2 eV thus at {two < leV) one is in low energy 
intraband response regime). The conductivity obtained 
is shown in Fig. (8) of their paper[l[ (it is also shown 
schematically in Fig. Hla) in the present manuscript). It 
is well known that spectrum below 1 eV in the cuprates 
cannot be fitted with single Drude model. This is also 
the case with CNCOC. It is this part of the spectrum 
that we will consider in our study. In their experi- 
mental study [H, they analyzed the experimental results 
with both two-component Drude model and generalized 
or extended Drude model (also called memory function 
formalismpBI]). We consider here their extended Drude 
model analysis of the experimental data. In this model 
all the spectrum below 1 eV is assigned to an itinerant 
state, with scattering rate and effective mass of charge 
carriers having frequency dependence. In this model the 
complex optical conductivity can be written as: 



ct(w) 



Aime 2 



47T \/t(lo) — iujm*(uj) 



(7) 



Here the scattering rate (T(uj) = T^jy) and effective mass 
m*(uj) both have frequency dependence. By fitting this 
with their experimental results they deduced frequency 
dependence of T(u;). They found that at low frequencies 
hit) < 0.4 eV scattering rate is almost proportional to 
ui (F(oj) = Tq + Cu). C turns out to be almost tem- 
perature independent and To increases with increasing 
temperature (ref. figure (10) in their paper[l|). They 
found that above hu> > 0.4 eV the scattering rate sat- 
urates to a constant value (very weakly dependent on 
temperature and doping). This is also shown schemat- 
ically in Fig. |U[b). They remark that this saturation 
behavior is similar to resistivity saturation. In the next 
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FIG. 4: Schematic view of the experimental results (see for 
original [J) 



section we will see that although the optical conductivity 
as calculated numerically using YRZ theory below 1 eV 
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appears qualitatively (to an unaided eye) in agreement 
with the experimental results (including the magnitude of 
conductivity (few hundred Vl~ 1 cm~ 1 )) but the bump 
in conductivity due to the presence of pseudogap 
shifts to lower energy regime with increasing dop- 
ing Ap g (x) = 0.6£o(l ~ x/0.2) in the numerical calcu- 
lation using YRZ theory and this is also reflected 
in the scattering rates. No doping independent 
saturation in T(lj) is observed in the numerical 
study. But in the experiments this kind of dop- 
ing independent saturation cut-off at haj ~ 0.4 eV 
in r(w) is observed. These are the main results of the 
present investigation and this is studied in greater detail 
in the next two sections. 



III. OPTICAL CONDUCTIVITY FROM YRZ 
ANSATZ AND COMPARISON WITH 
EXPERIMENT 

Our aim in this section is to numerically compute op- 
tical conductivity using YRZ ansatz for CNCOC sys- 
tem and to check how does it compare with the experi- 
ment. For the computation we needed a number of pa- 
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FIG. 5: (a) Optical conductivity (in Or 1 cm' 1 ) vs frequency 
(in eV) at x = 0.06 and at temperature T — 5K. Comparison 
at different 7 but with psuedogap. The value of the Drude 
scattering rate is given in the legend. All the energy scales are 
measured in eV . Figures in the second and 3rd row are for 
x = 0.08 and x — 0.1 respectively. In all figures to = 0.45eV . 



rameters of CNCOC system. These are taken from the 
literaturefl [ll, HI and tabulated below (Table I). 



TABLE I: Various parameters for CNCOC system. 



CitO-bond a- 
axis 


b-axis 


c-axis 


to 


~ 3.87 A 


~ 3.87 A 


~ 15 A 


~ 0.45 eV 


Hp at x = 0.06 


fj, p at x — 0.08 


a p at x — 0.1 




~ -0.2 eV 


~ -0.23 eV 


~ -0.25 eV 





The superexchange constant is J ~ 0.13 eV for 
Sr-iCuOzClz and for most of the cuprates to is in 
0.3 eV < t < 0.5 eV. Here we take t = 0.45 eV\v^. 
At x = 0.05, A° fl = 180 meV, ARPES do detect this 
pseudogap (Ap ! j iPES ~ 200 meV) at this doping around 
the point (wJ 0) in the brillion zone and also co-existing 
nodal metal 

For numerical computation we considered a 3-D sam- 
ple of 100 x 100 x 100 lattice points (thus with sample 
length, width, and hight: 100a, 100a, and 100c (where 
a is the a-axis (a = b) and c is the c-axis lattice con- 
stant)). This size of the sample is sufficiently large (as 
has been numerically verified) and size effects can be ne- 
glected. After performing the sum over k z in equation 
(3) we obtain 



<t(w,T) 



4ire 2 H 
N 2 ca 2 



+ I 7- + 1 , 



dyliy} L^L + ^ A(k x ,k y ,y)A(k x , k y ,u + y). 



Here we have redefined the units, now the frequency y, u) 
is measured in eV and {k x ,k y } = — -,— ^ + 6,...,+—, 
with S = N = 100. And the spectral function is 



A(k x , k y ,Lo) 



7 



{u> - e(k x ,k y ) - E pg (k x ,ky,u))) 2 + 7 2 



(9) 

Here 7 and lj measured in energy units (in eV). Band 
dispersion and self energy is also measured in eV and 
conductivity has the right units £l~ l cm~ 1 . 

Numerical computation is done on Mathematica-8 nu- 
merically by writing a small program. 

In Fig. ([5{a)) we plot optical conductivity as a func- 
tion of frequency at two different 7's and at temperature 
T = 5 K and hole doping x = 0.06. Dip seen at around 
0.2 eV is a signature of pseudogap which broadens with 
increasing 7 (red solid squares for 7 = 0.3 eV and blue 
filled circles for 7 = O.lto eV). This pseudogap signature 
(Dip at around 0.2 eV) disappears at higher tempera- 
ture. This is shown in Fig. @a)) where T = 300 K and 
we do not have pseudogap at this temperature. The con- 
ductivity obtained at 7 = 0.3 eV resembles closely what 
has been experimentally found (Fig. (8) of Waku etal.f^]). 
We will see in the next section that although it appears 
qualitatively in agreement with what has been found ex- 



6 



pcrimcntally but careful examination shows inconsisten- 
cies with YRZ. The other figures are plotted for x = 0.08 
and x — 0.1 as written in the figure titles. 



not seen in the experiment of Waku etal[l[. This is the 
subject of the next section. 
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FIG. 6: (a) With and without pseudogap at x = 0.06 and 
7 = 0.2 eV. Other plots are at different values of x and 7 as 
written in plot titles. 
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FIG. 7: With increasing doping x, size of the pseudogap de- 
creses Ap g (x) = 0.6io(l — x/0.2). This should result in left- 
ward shifting of the gap signature in the optical conductivity. 
Figure (a) show this shifting for x 1 — 0.06 and X2 = 0.08. 
Shift in A° g {x) is 6A° pg (x) = 3t (x 2 - xi) ~ 0.03eV. This 
roughly agrees with shift in the graph (the distance between 
x-axis ticks is 0.05eV). (b) Shifting for x = 0.06 and x = 0.1. 
Here also 6Ap„(x) = 3t (x 2 — Xi) ~ 0.054. This also roughly 
agrees with shift in the graph. 

In Fig. ([7(a)) we plot optical conductivity as a function 
of frequency at two different doping levels (x = 0.06, 
and x — 0.08) and at temperature T = 5 K and at 
7 = O.lto- There is a clear shifting of the pseudogap 
bump to lower frequencies with increasing doping x. This 
is consistent with decreasing A® vg (x) = 0.6io(l — x/0.2) 
with increasing x. In Fig. ((7jb)) the same shifting is 
shown for x = 0.06 and x = 0.1. Shifting seen in the 
graph agrees with SAp g (x) = 3t (x2 — x\) ~ 0.054 as 
explained in figure caption. This shifting of pseudogap 
signature in conductivity results similar shifting of the 
signature of pseudogap in the generalized Drude model 
scattering rate. But, in contrast, this kind of shifting is 



IV. EXTENDED DRUDE MODEL ANALYSIS 
AND INCONSISTENCIES WITH YRZ MODEL 



In this section we will show that optical conductivity 
as computed using YRZ is inconsistent with what has 
been experimentally observed inQ. We will deduce the 
above statement by the method of reductio ad absurdum 
of logic, i.e., proof by contradiction. So let us assume 
that YRZ is the correct model for the computation of op- 
tical conductivity in the low frequency regime which we 
are considering, i.e., the underlying transport properties 
of quasi-particles are captured by YRZ. Now, as done in 
the experimental paper by Waku etaif^J we analyze the 
optical conductivity from YRZ ansatz with the extended 
Drude model and extract the frequency dependent scat- 
tering rate. If the scattering rate so deduced agrees with 
scattering rate deduced with similar analysis of the exper- 
imental data, then, YRZ ansatz is consistent with what 
has been experimentally observed, otherwise, it is not. 

Now, as mentioned before (below equation (O), they[l[ 
deduce frequency dependence of the scattering rate T(co) 
of extended Drude model by using their experimental 
data. They found that at low frequencies hjj < 0.4 eV 
scattering rate is almost proportional to w (T(o;) = 
Tq + Ccu) and C turns out to be almost temperature 
independent. And To increases with increasing temper- 
ature (refer to figure (10) in their paper[l|). They see 
that above hiu > 0.4 eV the scattering rate saturates to 
a constant value. 

In our case, the scattering rate can be computed 
frorndl: 



r sd (w) 



r(w) 

9 

ne 



: Re \ 

\0~[UJ) 



cr(w) 



m* (cr(w)) 2 + (a lm (w)) 5 



(10) 



For this we need to compute imaginary part of conduc- 
tivity o~i m (uj) (real part is o~(uj)). As we are dealing with 
the casual perturbation and response relationship, the 
real and imaginary parts are related by Kramcrs-Kronig 
relations. By doing a numerical Kramcrs-Kronig inver- 
sion we obtain ai m {uj) (as plotted in Fig. ((Ha))) and 
then we calculate T g d- In Fig. (HJb)) we plot T(lo) = 

^vtL(^ (!VH = ^ r H) at 7 = 0.2 eV. 

We see clearly from Fig. (HJb)) that the scattering 
rate T(uj) obtained from extended Drude model mono- 
tonically decrease with increasing frequency. This is 
in sharp contrast with the experimental observations of 
Waku et al[l| where T(u>) first increases linearly with 
frequency and then saturates at around u> c ~ 0.4elA 
Also with reduced doping T(oj) shifts to lower hw in nu- 
merical study. Comparison of Fig. 2th) (or figure (10) 
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1 =0.2 eV, T = 5K 




(a) (b) 

FIG. 8: (a) Imaginary part of conductivity eri m (w) calculated 
from the real part by Kramers-Kronig transformation. The 
scattering rate 7 and temperature T is fixed at 0.2eV and 5K 
in these two graphs. The calculation is done for hole dop- 
ings x = 0.06, x = 0.08, and x = 0.1. We see that higher 
the doping higher is the conductivity and the peak in con- 
ductivity shifts to lower frequencies with increasing doping 
as is investigated in figure (|7j). The scattering rate T(u) ob- 
tained from extended Drude model monotonically decrease 
with increasing frequency. This is in sharp contrast with the 
experimental observations of Waku et al[l[ where T(lo) first 
increases linearly with frequency and then saturates at around 
uj c ~ 0.4eV. 




FIG. 9: The similar analysis as in figure 8 with 7 = 0.1. 



in[l|) with Fig. H^b) clearly points out the inconsisten- 
cies. The smeared out features of the pseudogap at 
higher 7 (7 = 0.2 eV in Fig. (jHJ) can be seen at lower 
7 (7 = 0.1 eV in Fig. ©). But the pseudogap peak in 
r(w) shifts to higher Hui with reduced doping (similar to 
the shift of ~ 0.05 eV between the peaks at x = 0.06 
and x = 0.1 in Fig. [3(b)). Also the behaviour of two 
scattering rates ( Fig.SJb) and Fig.^b)) is qualitatively 
different. 



V. CONCLUSION 



We theoretically investigated charge dynamics 
in weakly coupled Cu0 2 planes of the cupratc 
Ca 2 - x Na x Cu0 2 Cl 2 (CNCOC) using YRZ ansatz 
for the single particle Green's function in the pseudogap 
state. To an unaided eye, the results of our numerical 
calculation appears to be in qualitative agreement with 
what has been found experimentally in[l[. But a careful 
examination with extended Drude formalism shows that 
YRZ ansatz for the calculation of optical conductivity is 
not sufficient to understand the optical conductivity in 
Cu0 2 planes of the compound CNCOC. It seems that 
more physics is needed to fully understand the optical 
response especially the response from bound charges. 
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